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The momentum distribution of electrons from decays of heavy flavor (charm and beauty) for 
midrapidity \y\ < 0.35 in p + p collisions at y/s = 200 GeV has been measured by the PHENIX 
experiment at the Relativistic Heavy Ion Collider (RHIC) over the transverse momentum range 
0.3 < pt < 9 GeV/c. Two independent methods have been used to determine the heavy flavor 
yields, and the results are in good agreement with each other. A fixed-order-plus-next-to-leading- 
log pQCD calculation agrees with the data within the theoretical and experimental uncertainties, 
with the data/theory ratio of 1.72 ± 0.02 stat ± 0.19 sys for 0.3 < p T < 9 GeV/c. The total charm 
production cross section at this energy has also been deduced to be <r C c = 567 ± 57 stat ± 224 sys /ib. 

PACS numbers: 25.75.Dw 



Heavy-flavor (charm and beauty) production serves as 
a testing ground of QCD. Because of the large quark 
mass, it is expected that next-to-leading order pertur- 
bative QCD (NLO pQCD) can describe the production 
cross section of charm and beauty at high energy, partic- 
ularly at high pt- At the Tevatron, beauty production 
is well described by NLO pQCD 0. Charm production 
cross sections at high pt are found to be higher than the 
theory by ~ 50 %, but are compatible within the the- 
oretical uncertainties 0- Since heavy-flavor production 
at RHIC energies is dominated by gluon-gluon fusion, its 
production in polarized p + p collisions probes the gluon 
distribution G(x) and the gluon polarization AG(x). A 
good understanding of the reaction mechanism for heavy- 
flavor production is crucial for reliably extracting these 
distributions. Furthermore in Au+Au collisions at RHIC 
strong suppression of single electrons from heavy-flavor 
decays has been observed |3J Measurements of heavy- 
flavor production in p + p collisions provide a baseline for 
studying hot and dense matter effects in heavy ion reac- 
tions. Earlier measurements at RHIC 0,0 have a limited 
Pt range with substantial experimental uncertainties, so 
an improved measurement is crucial. 

We report the production cross section of electrons, 
(e++e~)/2, at mid-rapidity in p+p collisions at a/s = 200 
GeV for 0.3 < p T < 9 GeV/c measured by the PHENIX 
experiment. Contributions from semi-leptonic decays of 
heavy-flavor are determined using two independent meth- 
ods. This measurement has over two orders of magnitude 
larger statistics with much reduced systematic uncertain- 
ties compared to our previous measurement 

The data were collected by the PHENIX detector |fj| 
during the 2005 RHIC run using the two central arm 
spectrometers. Each spectrometer covers \i]\ < 0.35 
in pseudo-rapidity and A(f> — tt/2 in azimuth. It in- 
cludes a drift chamber (DC) and pad chambers (PCI) for 
charged particle tracking, a Ring Imaging Cerenkov de- 
tector (RICH) for electron identification, and an electro- 
magnetic calorimeter (EMCal) for electron identification 
and trigger. Beam-beam counters (BBCs), positioned at 
pseudo-rapidities 3.1 < \ri\ < 3.9, measure the position 
of the collision vertex along the beam (z v tx) and provide 
the interaction trigger. In this run, helium bags, one for 
each arm, were placed in the space between the beam 
pipe and DC to reduce multiple scattering and photon 
conversion. 



Two datasets are used for the analysis: (1) the min- 
imum bias (MB) dataset recorded by the BBC trigger, 
and (2) a "photon" trigger (PH) dataset triggered at 
level-1 requiring a minimum energy deposit of 1.4 GeV 
in an overlapping tile of 4 x 4 EMCal towers in coin- 
cidence with the BBC trigger. The PH trigger has ~ 
100 % efficiency for high pt electrons above 2 GeV/c in 
the active trigger tiles. The BBC trigger cross section is 
23.0 ± 2.2 mb. Since only ~ 50 % of inelastic p + p col- 
lisions satisfy the BBC trigger condition, only a fraction 
of the inclusive electron production events is triggered. 
This pt and process independent fraction is determined 
to be Cbias = 0.79 ± 0.02 from the yield ratio of high pt 
w° with and without the BBC trigger. After selection of 
good runs and a vertex cut of |z vtx | < 20 cm, an inte- 
grated luminosity (£) of 45 nb _1 in the MB dataset and 
1.57 pb _1 in the PH dataset are used for the analysis. 

Charged particle tracks are reconstructed using DC 
and PCI and confirmed by a hit in EMCal within 4cr 
in position. The momentum resolution is cr p /p ~ 0.7% © 
0.9%p (GeV/c), and the momentum scale is calibrated 
within 1% using the reconstructed mass of J/ip — > e + e~ . 

Electron identification (elD) requires at least two as- 
sociated hits in RICH, a shower shape cut in EMCal, 
and a cut in the ratio E/p where E is energy mea- 
sured in EMCal. We require 0.7 < E/p < 1.3 for 
0.8 < pt < 5 GeV/c. For lower pt, the minimum 
value of E/p decreases with decreasing pt to 0.55 at 
Pt = 0.3 GeV/c. The E/p cut removes background elec- 
trons from photon conversions and semi-leptonic decay 
of kaons (K — > evir(K e z)) that occur far from the ver- 
tex, and most of the remaining hadron background. The 
hadron contamination after the E/p cut is 3% at pt = 0.3 
GeV/c and less than 1% for 0.8 < p T < 5 GeV/c with 
elD efficiency of approximately 90 %. 

For pt > 5 GeV/c, where pions also emit Cerenkov 
photons in RICH, tighter electron identification cuts are 
applied. We require at least 5 associated hits in RICH, a 
tighter shower shape cut in EMCal, and 0.8 < E/p < 1.3. 
With these cuts, the electron measurement is extended 
to 9 GeV/c in pt- The elD efficiency of the tighter cuts 
is pt independent, and is determined to be 57% of that 
for px < 5 GeV/c by applying the same tighter cuts for 
Pt < 5 GeV/c. With the tighter cuts, hadron contami- 
nation is negligible for pt < 7 GeV/c. For 7 < pt < 8 
(8 < p T < 9) GeV/c, a 20% (40%) hadron contamina- 
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tion is determined and subtracted using a Gaussian plus 
exponential fit to E/p distribution. 

The invariant cross section for electron production is 
calculated using the following formula, 

d 3 a 1 1 N e 1 f 
dp 3 C 2np T Ap T Ay e roc e b i as ' 

where N e is the measured electron yield; e rcc , calculated 
using a full GEANT 7] simulation, includes the geomet- 
rical acceptance, track reconstruction and cID efficiency, 
and the smearing effect due to finite momentum reso- 
lution. For the PH dataset, e rcc also includes the PH 
trigger efficiency. The cross sections from the MB and 
the PH datasets are consistent with each other for the 
overlapped pt region. 

The inclusive electron yield consists of three compo- 
nents: (1) electrons from heavy-flavor decay, (2) "pho- 
tonic" background electrons from Dalitz decays of light 
mesons and photon conversions primarily in the beam 
pipe, and (3) "non-photonic" background electrons from 
the remaining K e % decays and dielectron decays of vector 
mesons. The photonic background is much larger than 
the non-photonic background. We determined the spec- 
trum of electrons from heavy-flavor decay by subtracting 
the background components from the inclusive spectrum 
using the following two independent methods. 

In the "cocktail subtraction" method 0, 0- El a cock- 
tail of electron spectra from various background sources 
is calculated using a Monte Carlo event generator of 
hadron decays. The most important background is the 
7T° Dalitz decay, so we use our measured 7r° and tt^ spec- 
tra as input to the generator. The spectral shapes of 
other light hadrons h are obtained from the pion spec- 
tra by hit scaling (pt — > \/.Pt + M£ — M%). Within 
this approach the ratios h/ir° are constant at high pt- 
For the relative normalization, we use the following ra- 
tios: rj/ir = 0.48 ± 0.03 9], P°/tt° = 1.0 ± 0.3, w/vr° = 
0.90 ±0.06 [13, rf/TV = 0.40 ±0.12, </>/tt° = 0.25 ±0.08. 
For pt > 2 GeV/c, contributions from r\ and all other 
hadrons combined are approximately 20% and 10% of ir°, 
respectively. Another major background electron source 
is conversions of photons in the beam pipe (0.29% of a 
radiation length(Xo)) as well as in the air and the he- 
lium bags (0.1% Xo). The conversion electron spectrum 
is very similar to that of Dalitz decays. Using a detailed 
GEANT simulation of the PHENIX detector, the ratio 
of electrons from conversions to Dalitz decays, i?cD, is 
determined to be 0.40 ± 0.04 for 7r°, essentially gj inde- 
pendent. i?cD is approximately half of that in jj| since 
the helium bags eliminated most of the conversions out- 
side of the beam pipe. The conversion spectra are calcu- 
lated by scaling the Dalitz decay spectra by Rcd, with 
small corrections to account for the species dependence 
of the relative branching ratio of Dalitz decay to photon 
decay ((h — > ee^)/(h — > 77)). The internal and exter- 
nal conversions of direct photons are also included in the 



cocktail, using our measured direct photon spectrum |ll| 
as input. The direct photon contribution is comparable 
to or greater than that from the ij for pt > 5 GeV/c. 
Non-photonic backgrounds are also included in the cock- 
tail. Since the K e ^ background depends on the analysis 
cuts, it is evaluated by a full GEANT simulation. 

In the "converter subtraction" method ^^1, we in- 
troduce an additional photon converter (a thin brass 
sheet of 1.67% Xq) around the beam pipe for part of 
the run. The converter multiplies the photonic electron 
background by a fixed factor, i? 7 ~ 2.3, which is de- 
termined precisely via GEANT simulation. i? 7 is larger 
than in [12| since we have less conversion material in the 
2005 run. The photonic background N2 is determined as 
AH = (A e c -(l-e)A e NC )/(i? 7 -l+e), where A e c and A e NC 
are electron yield with and without the converter, respec- 
tively; and e (2.1%) represents a small loss of electrons 
due to the converter. The non-photonic component is 
then determined as A e non -T = A e NC - N]. Small remain- 
ing non-photonic background, such as K e ^ and hadron 
contamination, are subtracted. 

These two methods are complementary to each other. 
The converter method is more accurate, and it allows us 
to extract a heavy-flavor signal down to pt = 0.3 GeV/c 
where the signal is only w 10% of inclusive electrons. 
In addition, the measured photonic background is 
used to confirm and to calibrate the normalization of the 
calculated cocktail yields. A drawback of the method is 
statistical precision: the converter run contains only a 
small fraction (~ 7% in the 2005 run) of the data. The 
cocktail method can use the full statistics at high pt, 
where the photonic background becomes a small fraction 
of inclusive electrons. 

Systematic uncertainties are categorized into (a) in- 
clusive electron spectra, (b) cocktail subtraction, and 
(c) converter subtraction. Category (a) is common to 
both analyses, and includes the uncertainties in luminos- 
ity (9.6%), geometrical acceptance (4%), elD efficiency 
(3%), and the PH trigger efficiency (3% at the plateau). 
Uncertainties in cocktail subtraction (category (b)) in- 
clude the normalization (8%) and pt dependent shape 
uncertainty (2% at pt — 2 GeV/c, increasing to 6% at 
9 GeV/c). In the converter analysis (category(c)) the 
dominant uncertainties are in i? 7 (2.7%) and in the rel- 
ative acceptance in the converter and the normal runs 
(1.0%). These uncertainties are propagated into the un- 
certainties in the heavy-flavor electron yields and added 
in quadrature. 

Figure shows the ratio of the measured to the 
cocktail calculation as a function of pt- The ratio is con- 
sistent with unity within the uncertainties of the cocktail. 
At high p T (> 1.8 GeV/c), the ratio is 0.94 ± 0.02 stat on 
average. Since this is within the uncertainty of the cock- 
tail normalization, we rescale the cocktail yields by this 
factor. This removes the 8% normalization uncertainty 
in the cocktail. 
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FIG. 1: Ratio of photonic electrons measured by the con- 
verter method to the cocktail calculation. Data from the MB 
(PH) dataset are shown below (above) 1.8 GeV/c. The upper 
and lower curves show the systematic error of the cocktail. 
Error bars are statistical only. 
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FIG. 2: Ratio of non-photonic electrons to photonic back- 
ground. Error bars are statistical errors and the error bands 
show the cocktail systematic errors. The solid, dashed, dot- 
dashed, and dot curve is the remaining non-photonic back- 
ground from K e s, p — > ee, uj — > ee, and hadron contamina- 
tion, respectively. 



In Fig. filled circles (squares) show the ratio of non- 
photonic electrons relative to photonic background de- 
termined by the converter (cocktail) method. The non- 
photonic electrons are dominantly heavy-flavor decay sig- 
nals. The remaining non- photonic background contribu- 
tions have been calculated and are shown in Fig. |2J The 
two methods are consistent with each other. The ra- 
tio monotonically increases with increasing px, becoming 
greater than unity for px > 2.4 GeV/c, and saturates at 
~ 3 for pt > 5 GeV/c. The large signal-to-background 
ratio is due to the small amount of conversion material 
in the spectrometer acceptance. 

Figure [3] (a) shows the invariant differential cross sec- 
tion of electrons from heavy-flavor decays. All back- 
ground has been subtracted, including the non-photonic 
components shown in Fig. [3 The data from the two anal- 
ysis methods are combined: at lowpx ipr < 1-6 GeV/c) 
the converter subtraction method on the MB data set is 
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FIG. 3: (a) Invariant differential cross sections of electrons 
from heavy-flavor decays. The error bars (bands) represent 
the statistical (systematic) errors. The curves are the FONLL 
calculations (see text), (b) Ratio of the data and the FONLL 
calculation. The upper (lower) curve shows the theoretical 
upper (lower) limit of the FONLL calculation. In both panels 
a 10% normalization uncertainty is not shown. 



used; at intermediate pt (1.6 < pt < 2.6 GeV/c) the 
converter method on the PH dataset is used; and at high 
Pt (pt > 2.6 GeV/c) the cocktail method on the PH 
dataset is used. 

The data are compared with a fixed-order-plus-next- 
to-leading-log (FONLL) pQCD calculation HaHH. The 
top curve in Fig.|3shows the central values of the FONLL 
calculation. The contributions of charm and beauty are 
also shown. For px > 4 GeV/c, the beauty contribu- 
tion becomes dominant. In Fig. [3] (b), the ratio of the 
data to the FONLL calculation is shown. The ratio is 
nearly px independent over the entire pt range. Fitting 
to a constant for 0.3 < pt < 9.0 GeV/c yields a ratio of 
1.72 ± 0.02 stat ± 0.19 sys Similar ratios are observed in 
charm production at high px at the Tevatron [2| . The up- 
per limit of the FONLL calculation is compatible with the 
data. Recently STAR reported ^| that non-photonic 
electron production in p + p at y/s = 200 GeV is 5.5 
times larger than predicted by the same FONLL calcu- 
lation. We do not observe such a large discrepancy. We 
note that the photonic electron background in our spec- 
trometer is approximately 1/10 that of STAR due to the 
small amount of conversion material in the PHENIX ac- 
ceptance. 

The total charm cross section is derived by integrat- 
ing the heavy-flavor electron cross section for px > 0.4 
GeV/c: da e (p T > 0A)/dy = 5.95 ± 0.59 ± 2.0 fjb. 
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The systematic error is obtained by integrating the up- 
per and lower systematic error limits of the differen- 
tial cross sections, since the systematic errors are es- 
sentially coherent. The cross section is then extrapo- 
lated to px = using the spectrum shape predicted 
by FONLL: da e (p T > 0)/dy = 10.9 ± 1.1 ± 3.8 fib. 
We have assigned 10% to the systematic uncertainty 
of the extrapolation, and have subtracted contribution 
from beauty and beauty cascade decays (0.1 fib). We de- 
termine the charm production cross section, da cS /dy = 
123 ± 12 ± 45 fib, by using a c — > e total branching ra- 
tio of 9.5 ± 1.0%, calculated using the following charmed 
hadron ratios: D+ /D° = 0.45 ±0.1, D s /D a = 0.25 ±0.1, 
and A c /L>° = 0.1±0.05. The rapidity distribution of elec- 
trons is broader than that of D mesons due to decay kine- 
matics. A correction to this effect (7%) has been applied. 
Using the rapidity distribution from HVQMNR [l(j with 
CTEQ5M 17] PDF, the total charm cross section is de- 
termined to be cr cE = 567 ± 57 stat ± 224 s y s fib. We have 
assigned 15% systematic error to the extrapolation. This 
result is compatible with our previous measurement 
(920 ± 150 ± 540 fib) and the value derived from Au+Au 
collisions H3 (622±57±160 fib per NN collisions). The 
FONLL cross section (256i^g fib) is compatible with the 
data within its uncertainty. STAR has reported a some- 
what larger value in d+Au 5] (1.3 ± 0.2 ± 0.4 mb 
per NN collisions). Although the data extend to high 
Pt where the beauty contribution is expected to be dom- 
inant, the present analysis does not separate charm and 
beauty contributions. The beauty cross section predicted 
by FONLL is 1.87±g;|? fib, and the upper FONLL curve 
is consistent with the data. 

In conclusion, we have measured single electrons from 
heavy- flavor decays in p + p collisions at y/s — 200 GeV. 
The new data reported here provide a crucial baseline for 
the study of heavy quark production in hot and dense 
matter created in Au±Au collisions. The agreement be- 



tween the data and the FONLL pQCD calculation within 
the theoretical and the experimental uncertainties sug- 
gests that a reliable extraction of gluon polarization from 
heavy- flavor production in polarized p+p collisions is at- 
tainable. 
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